Mon. Not. R. Astron. Soc. 000, [HI] (2010) Printed 19 June 2012 (MN WT^ style file v2.2) 



Variable Accretion Rates and Fluffy First Stars 



O 

00 



Rowan J. Smith^*, Takashi Hosokawa^, Kazuyuki Omukai^, 
Simon CO. Glover \ and Ralf S. Klessen^ 

^ Zentrum fiir Astronomie der Universitdt Heidelberg, ITA, Albert-Ueberle-Str. 2, 69120 Heidelberg, Germany 
^ Jet Propulsion Laboratory, California Institute of Technology, Pasadena CA 91109, USA 
^ Department of Physics, Kyoto University, Kyoto 606-8502, Japan 



19 June 2012 



o 

o 

in: 

(N 
> 

in 



ABSTRACT 

We combine the output of hydrodynamical simulations of Population III star cluster 
formation with stellar evolution models, and calculate the evolution of protostars 
experiencing variable mass accretion rates due to interactions within a massive disk. 
We find that the primordial protostars are extended 'fluffy' objects for the bulk of 
their pre-main-sequence lifetimes. Accretion luminosity feedback from such objects is 
high, but as shown in previous work, has a minimal effect on the star cluster. The 
extended radii of the protostars, combined with the observation of close encounters in 
the simulations, suggests that mergers will occur in such systems. Furthermore, mass 
transfer between close protostellar binaries with extended radii could lead to massive 
tight binaries, which are a possible progenitor of gamma ray bursts. 
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1 INTRODUCTION 



It is becoming increasingly clear that Population III star for- 
mation is not as different from that in the present day as was 
once thought. Several authors have now shown that instead 
of a single massive star being formed at the centre of a pri- 



mord i al halo, a sm all cluster is formed ( e.g. iMachida et al 



200S: IClark et al 



2010l : iPrieto et al 



2011a,b; IXurk et all |2009| : IStacv et al 



I2OIII ). There are frequent dynamical in- 



teractions within these clusters that can strongly alter the 
accretion rates of the protostars l|Greif et al.ll201ll ). As such. 
Population III star formation begins to resemble simulations 
of local clustered star formation (e.g Bonnell et al.l I2OIII : 
lOffner et all I2OI0I : iGirichidis et al.l |201l') . For a review of 
presen t day clustered star formation see Mac Low fc KlessenI 
l|2004 ). 

A parallel strand in the study of Population III star 
formation has been the investigation of the stellar evolu- 
tion of the protostar as it grows in mass. As the protostar 
grows, the consequences of its protostellar evolution become 
increasingly significant. For example, during the protostel- 
lar stage of the evolution, radiation from the protostar alters 
the thermal and chemical structure of the circumstellar disk. 
As a result, the nature of disk fragmentation via the grav- 
itational inst ability is affected by the p rotostellar radiation 
l|Clark et al.l 2011b; ISmith et adbOUl ). As the star grows 
in mass and approaches the main sequence, it also starts to 
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emit UV radiation. This creates an HII region which can pre- 
vent further mass accretion onto the star and limit t he final 
stellar mass (iMcKee fc TanI l2008l : iHosokawa et all 2011b; 
IStacv et al.ll2011 ). The stellar evolution also determines the 
protostellar radius, which is an important quantity for un- 
derstanding the dynamic evolution of a young cluster, as 
this determines the likelihood of protostellar collisions and 
me rgers (for an analysis of pres ent day protostellar clusters 
see iBaumgardt fc Klessenir201ll ). 

Much of the pioneering work on the stellar evolu- 
tion o f Popu lation III protost ar s was done b y IStahler et al.l 
l|l986l) and lOmukai fc Pallal (1200 ll . 120031 ) who numeri- 
cally solved the interior structure equations assuming var- 
ious constant accretion rates and showed that the evo- 
lution quantitatively changes with different rates . How- 
ever, as was pointed out in ISmith et al. I (|201ll . here- 
after Sll), the protostars in a Population III minihalo 
do not accrete at a constant rate, but at a highly vari- 
able one, as the environment from which they accrete is 
not uniform. Time-variable mass accre tion is also ubiqui- 
tous i n present-day star f ormat i on (e.g. Schmeia fc KlessenI 



|2004 IVorobvov fc Basul l2006l : iMachida et all |2011^ and 

its effects on protostellar evolut i on h av e begun to 
be e xplored ("Wuchterl 
20091: iBaraffe fc Chabrier 



Klesseh '200ll ; [Saxa-ffe et all 
2010; Hosok awa et all 2011a; 



Hartmann et al.l l201ll) . In this paper, we seek to combine 



the strands of hydrodynamic simulations of fragmentation 
with stellar modelling and study how time- variable mass ac- 
cretion affects the evolution of primordial protostars. 
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2 PROTOSTELLAR ACCRETION RATES 

We use protostellar accretion data taken from the sim- 
ulations of Sll. These are smoothed particle hydrody- 
namics (SPH) r e-simulations of f ive cosmological miniha- 
los presented in iGreif et al.l (|201ll ') in which the heating of 
the gas by the accr etion luminosity produced by accret- 
ing 'sink particles' jSate et alj Il995l : IJappsen et al.l l2005l : 
iFederrath et al.ll20ld ) has been taken into account. Sink par- 
ticles represent regions where protostars will form. They are 
non-gaseous particles that interact with their environment 
only through gravitational forces and are formed from bound 
and unambiguously collapsing fragments which are above 
a specified critical density, often chosen to be the density 
at which the Jeans length of the gas is only just resolved. 
Through the representation of protostars as sink particles, 
the mass going into star formation, and the location of pro- 
tostars can be tracked over the lifetime of the simulation. 
Relatively large sink radii of 20 AU were used in Sll, which 
means that the inner disks surrounding our protostars are 
not resolved. However, by sacrificing the ability to follow 
these regions the computational time for the calculation is 
reduced sufficiently to allow the evolution of the protostars 
to be followed up to the point where they start contracting 
to the main sequence and p roduce substantial ionising ra- 
diation (|Tan fc McKeel I2OO8I ) . Our large sink particles also 
slightly suppress the total number of stars formed in our 
simulation, as we can no longer form fragments in the inner 
disk. These simulations, by having large sinks and by in- 
cluding accretion luminosity feedback, therefore represent a 
conservative estimate of the total amount of fragmentation 
that will occur within primordial halos. 

In Figure [1] we show two snapshots of the central region 
of Halo 4 from Sll. An extended thick disk is formed at the 
centre of the halo which h as strong spi ral arm features due 
to gravitational instability. Clar k et al.l (2011b) showed that 
primordial accretion disks become unstable due to the rate 
at which material can be transported inwards through the 
disk being lower than the rate at which new material falls on 
to the disk. This effect causes the disk in Figure [T] to frag- 
ment into multiple additional sinks which orbit around the 
centre of mass of the system within the irregular accretion 
disk. The calculations in Sll had to be terminated when 
protostars had masses in excess of 20 Mq as they did not 
include the effects of ionisation feedback. However in reality 
the protostars will continue to accrete beyond this point, 
until the ionising radiation finally quenches the accretion 
iHosokawa et al.l 2011b'). 

The accretion rate of the orbiting sink in Figure [1] is 
shown in the upper panel of Figure (2] The accretion rate is 
highly variable with sharp peaks and troughs. Figure[T]shows 
that the sink orbits around the protostar at the centre of 
the disk. When the protostar is at the periastron of its orbit 
(as shown in the top panel of Figure [T]) it passes though 
a dense one-armed spiral that causes the protostellar disk 
that surrounds the sink to grow in size, leading to a high 
accretion rate. On the other hand, when the protostar moves 
through the low density material outside the spiral arms (as 
shown in the bottom panel of Figure [l}, its protostellar disk 
becomes depleted of gas, and the sink has a lower accretion 
rate. This variation leads to fluctuations of up to an order 
of magnitude in the accretion rate. 
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Figure 1. A column density projection of the central regions of 
halo 4 from Sll where there is a gravitationally unstable disk from 
which multiple protostars are formed. Black circles show the loca- 
tion of sink particles representing protostars. Over-plotted on the 
figure is the path of one of the sink particles (shown in blue) which 
is orbiting within the disk. The central sink is shown in green. The 
top panel shows when the orbiting sink has a high accretion rate 
(M = 0.014 M0yr-l, M* = 10.35 Mq) as it passes through 
the spiral arm close to the central sink and the bottom panel 
shows where it has a low accretion rate (M = 0.002 MQyr~^, 
M» = 10.45 Mq) as it passes through a less dense region. 



In our simulations we are missing the inner 20AU of 
the accretion disk as it is inside the sink particle radius. This 
raises the possibility that some of the variability in our mea- 
sured accretion rates might have been smoothed by an inner 
accretion disk. However, as shown in IClark et all (2011b), 
primordial accretion disks are massively self-gravitating and 
rapidly accreting even down to scales of 1.5 AU. Further- 
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Figure 2. Top: The accretion rate of the orbiting sink in Fig- 
ure [T] The sohd red hne shows the accretion rate measured di- 
rectly from the simulation. The mean accretion rate at which 
2OM0was accreted (2.4 X 10~^ M0 yr"-*^) is shown by the dot- 
ted blue line and a running average of the variable rate is shown 
by the dashed green line. An additional higher mean accretion 
rate of 6 X 10~^ Mq yr~^ is plotted for comparison using a pink 
dot-dashed line. Bottom: The evolution of the protostellar radius 
given the above accretion histories. The variable rate measured 
from the simulation results in the protostar having an extended 
stellar radius which is sustained over a large mass range. 



more, the simulations in Sll did not use the instantaneous 
accretion rate, but instead used a smoothed average calcu- 
lated over the last 100 years. This should to some extent 
mimic the effect of accreted material being buffered by the 
inner disk. 

3 PROTOSTELLAR EVOLUTION WITH 
VARIABLE ACCRETION RATES 

To calculate the stellar radii we use the measured variable 
accretion rates as in put for a stellar evolution code, an ap- 
proach first used in IWuchterl fc KlessenI (|200ll ) . We calcu- 
late the interior structure of the protostar by numerically 
solving the four stellar structure equations (i.e. the equa- 
tions of continuity, hydrostatic equilibrium, energy conser- 
vation and energy transfer) with the inclusion of mass ac- 
cretion ( Omuka i fc P alla 2003; Hosokawa & Omukai 2009; 
iHosokawa et aLlboiol '). 

The thermal efficiency of mass accretion, i.e., the spe- 
cific entropy of the accreting material, is an unknown but 
important quantity which controls the effect of mass accre- 
tion on the structure of the protostar. Our numerical code 
handles the two opposite limits where the mass accretion 
is thermally efficient ( "hot accretion" ) and inefficient ( "cold 
accretion" ) by ad opting different outer boundary conditions 
l|Hosokawa et al.l 2011a). The exact value of the thermal ef- 



Figure 3. As Figure [2] but for the central sink. In this case 
the mean accretion rate is 1 X 10~^Mq yr~^. The high initial 
accretion rate of the protostar, before fragmentation limits its 
accretion, causes the protostellar radius to be enhanced compared 
to the averaged rate. 



ficiency of mass accretion should depend on the thermal 
structure of the accretion flow connecting the protostar and 
the disk. It is naively expected that the accretion flow be- 
comes hotter for a higher accretion rate. Studies of FU Ori- 
onis outbursts in the Galaxy, for which mass accretion at 
~ 10~* Mq yr~^ is expected, suggest that the gas tem- 
perature near the inn er edge of the disk is ~ 10* K (e.g., 
iHartmann et al.ll201l[ ). The typical accretion rates in Pop- 
ulation III star formation are 10"'' Mq yr~^ or more, and 
the disk temperature ne ar the stellar surface should exceed 
lO'^ K in such situations ()McKee fc Tanll2008l 'l. As such, hot 
accretion is the most realistic case for our model and we use 
this assumption in our following calculations. 

Using our stellar evolution code we calculate the proto- 
stellar radius using four rates: the measured accretion rate 
from Sll, a higher accretion rate, the mean accretion rate, 
and a running average of the measured rate. The mean ac- 
cretion rate is found by dividing the final mass of the pro- 
tostar in Figure [2] by the protostar's age. The running aver- 
age is calculated by averaging the measured accretion rate 
over a ±3 Mq mass range. This has the effect of smoothing 
out short period variations of < 10^ years in the accretion 
rate. Note that the simulations follow the evolution roughly 
over 10* years. The overall trend that the accretion rate 
decreases as the stellar mass increases is extracted in the 
running-average case. 

The resulting stellar radii are shown in the bottom panel 
of Figure (2] The variable accretion rate has increased the 
radius of the protostar compared to the mean rate. The 
stellar radius now extends to a maximum value of over 
200 Rq , or roughly lAU. The maximum value of the proto- 
stellar radius before it reaches the main sequence is similar 
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Table 1. The peak and mean radii for the measured variable 
accretion rates and the constant mean rate. 





Variable rate 


Mean rate 


Orbiting protostar 








Peak Radius 


255.7 


Rq 


166.7 R© 


Mean Radius 


102.5 


Rq 


63.4Rq 


Central protostar 








Peak Radius 


165.0 


Ro 


108.8 R© 


Mean Radius 


107.5 


R© 


35.9 R0 



to what we find when we use a higher constant accretion 
rate of 6 x 10"'^ M0 yr~^. This is unsurprising, as we see 
that the running average of the measured accretion rate is 
roughly equal to 6 x 10"'^ M© yr~^ for protostellar masses 
M* < 12 Mq. For larger M, , the running-average of the 
accretion rate decreases significantly, eventually becoming 
comparable to the mean accretion rate. Because of this de- 
crease, the use of the higher constant accretion rate does 
not yield an accurate prediction for the protostellar radius 
at late times. After the protostar starts contracting to the 
main sequence, the stellar radii of the mean, running-average 
and full models converge. Consequently, the end state of the 
protostar is set by the mean accretion rate and is not af- 
fected by the earlier variability in its accretion history. 

As the variation shown for the protostar in Figure [2] is 
an extreme case, we also repeat this calculation with the pro- 
tostar at the centre of the disk in Figure [l] Figure [3] shows 
that once again the protostellar radius is substantially en- 
hanced relative to the mean rate. The accretion rate of the 
central protostar drops rapidly after it reaches a mass of 
around 6Mq. This corresponds to the point in our simula- 
tion where the disk around the star becomes unstable and 
the additional protostars in the halo are formed. The new 
protostars also accrete gas from the massive central disk and 
hence reduce the accretion rate of the central ob ject in a pro- 
cess termed fragmentation-induced starvation fPet ers et al.l 
[2010 ). The evolution of this protostar's radius cannot be re- 
produced by any model with a constant accretion rate as it 
has two peaks in its evolution and a constant accretion rate 
can only produce one. 

Table [T] shows the maximum and average value of the 
stellar radius of each sink for the variable accretion rate and 
its mean. In both cases the protostellar radius has a higher 
maximum and a higher mean value when the variable rate 
is used. In particular, for the central sink shown in Figure [3] 
the mean protostellar radius is three times higher than that 
found using the mean accretion rate. 

We find that the young protostars have large initial radii 
due to having an enhanced rate of accretion during the early 
phase of their evolution relative to that in later stages. In 
Figure [2] the accretion rate of the sink is at an initially high 
level due to the sink accreting large amounts of material 
from the spiral arms within the dense disk. In Figure [3] the 
rate was initially high but the formation of protostars in a 
disk around the central protostar caused the accretion rate 
in the later stages to fall due to fragmentation-induced star- 
vation. Additional dynamical ejections can also remove sinks 
from the dense central regions of the halo and cause their 
accretion rates to fall, which once again means that the ac- 



cretion occurs disproportionately in the early stages of the 
protostars' evolution. 

The importance of the initial accretion rates can be 
seen by considering the phases of stellar evolution. While 
the protostellar mass is less than around 10 M©, the proto- 
stellar evolution is adiabatic, as the accretion timescale for 
the protostar (i.e., the timescale on which it increases its 
mass) is shorter than the Kelvin-Helmholtz (KH) timescale 
(i.e., the timescale on which the star can redistribute entropy 
internally) 



where M, is the protostar's mass, J?, the radius and L* 
the luminosity. Rapid accretion during this stage dumps 
large amounts of entropy onto the surface of the proto- 
star, which cannot be quickly redistributed, causing the stel- 
lar surface to swell. The opacity in the stellar interior de- 
creases as the stellar mass increases, which in turn cause s 
the stellar luminosity to increase (e.g., IStahler et al]|l986l ). 
The KH timescale then becomes shorter and shorter, and 
finally falls below the accretion timescale. After that point 
(M* ~ 12 Mq in both the cases), the protostar begins to 
contract, losing its energy via radiation. The characteris- 
tic timescale of the evolution is the accretion timescale in 
the early adiabatic accretion stage and the KH timescale in 
the later contraction stage. Both timescales are around 10^ 
years. 

Using the running average of the accretion rate allows 
us to investigate whether the short-term variability or the 
long-term trend in the accretion rate is the greater factor 
in determining the stellar evolution. In Figures [2] and [31 the 
protostellar radius calculated from the running average of 
the accretion rate provides an excellent match to the true 
value. Consequently, the long-term trend of the accretion is 
more important than the short-term variability. It is reason- 
able that accretion variability shorter than the characteristic 
timescale of ~ 10'^ years does not significantly influence the 
overall evolution of the protostar. 

Generally protostellar accre tion rates are predicted t o 
fall exponentially with time (e.g. ISchmeia fc Klessen|[2004l '). 
but for the orbiting protostar the accretion rate is main- 
tained at a high value throughout its early evolution which 
increases its protostellar radius. Conversely, the accretion 
rate experienced by the central protostar rapidly decreases 
after the disk fragments but during the adiabatic stage of its 
evolution the accretion rate is very high (over 0.01 Mq yr~^) 
and consequently it has a large initial stellar radius. The 
process of disk fragmentation therefore leads to a system in 
which the protostars have extended radii. 

iHosokawa fc Omukail (|2009l ) have shown that all mas- 
sive protostars accreting at high rates are fluffy extended 
objects, but the Population HI stars here are particularly 
striking examples. While the final main sequence radius of 
the stars considered here is small, the large radius during 
the protostellar stage has ramifications for the type of stellar 
cluster that is formed within the primordial minihalo. This 
is due to the accretion rate and stellar radius of the proto- 
star determining both the strength of protostellar feedback 
in the young cluster, and the likelihood of stellar interactions 
and mergers, as we shall discuss in the following sections. 
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4 FEEDBACK 

The dominant form of stellar feedback during the first few 
thousand years of a protostar's life in a primordial halo is 
that from accretion luminosity. The luminosity is given by 



Lace — 



GA'LM 
R, ' 



(2) 



where M is the accretion rate. The accretion rate contributes 
doubly to this relationship, as in addition to the direct de- 
pendence of Lace on M, the accretion rate also affects the 
radius of the star. Figure|4]shows the accretion luminosity of 
the previously considered orbiting sink calculated with the 
measured accretion rate, the constant rate and the decreas- 
ing accretion rate. During the early evolution of the proto- 
star the variable accretion rate gives accretion luminosities 
above 10'* Lq, which is substantially higher than the value 
we find for the constant accretion rate. Sll considered the 
effects of accretion luminosi ty in detail using a sem i-analytic 
model based on the work of lOmukai fc Palla (|2003l ) for radii 
and accretion luminosities. In this semi-analytic model the 
stellar radius evolves as 



R. 



26M,"2^(M/10-=^)°- 



.41 



M, ^ Pi 

Pi ^ M, < P2 

P2^M, &i R<R„ 



(3) 



where J?* has units of R©, M* has units of M©, and the 
accretion rate M has units of M0 yr~^. In this expression, 
the transition points between the evolutionary phases, pi 
and p2, scale as 



Pi = 5(M/10-3)°-"^ Mq, 
P2 = 7{M/W'^y'-^^ Mq, 

and the main sequence radius is given by 



(4) 



(5) 



The accretion luminosity is then calculated from the mea- 
sured accretion rate and the stellar radius. This model is 
plotted in Figure 3] in black, and is a good match to the true 
accretion luminosity, particularly during the early stages of 
the protostar's evolution where the star evolves adiabati- 
cally. We therefore find that the conclusion of Sll that ac- 
cretion luminosity does not prevent fragmentation of primor- 
dial gas is valid, despite the fact that Sll used a simplified 
semi-analytic approach rather than a full stellar evolution 

model. 

As shown in lHosokawa et al] (2011b), ionising radiation 
plays a key role in determining the final mass of a Population 
III protostar. We find that the protostar reaches the main 
sequence at the same stage with either a constant or vari- 
able accretion rate, and therefore in both cases the protostar 
starts ionising its surroundings at the same mass. 



5 CLOSE ENCOUNTERS 

As the protostars are formed from disk fragmentation, the 
separation between sink particles is necessarily small. In 
fact, close encounters seem to be a feature of the small 
clusters formed in Population III halos. In the study of 
iGreif et al.l l|201ll ). which considered fragmentation on scales 
down to 100 Rq, many encounters were seen where two sink 
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Figure 4. A comparison of the accretion luminosity from the 
previously considered orbiting sink. The solid red curve uses the 
variable accretion rate, the dotted blue curve the constant rate 
and the dashed green curve is the running average rate. The black 
dashed-dotted curve shows the accretion luminosity calculated 
using the semi-analytic method used in Sll. 
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Figure 5. The separation between the sink particles in the tight- 
est binary seen in our hydrodynamic simulations (black). The 
red line shows the total stellar radius of the protostars during the 
evolution. As the radii of the protostars is frequently compara- 
ble to their separation these particular sinks are no longer well 
represented by point masses in this extreme case. 



particles came within this radius. In the simulations consid- 
ered here, we do not resolve fragmentation at separations 
of less than 20 AU and so do not follow the evolution of 
fragments formed in the inner parts of protostellar accre- 
tion disks. We therefore see fewer of these close encounters. 
Despite this, in all but one of the five halos studied in Sll 
there is at least one pair of sink particles which come within 
a distance of 1 AU of each other. 

As an illustrative example, in Figure [5] we plot the dis- 
tance between two sinks that make up a tight binary in halo 
3 of Sll. This binary is the closest found in the simulations 
carried out by Sll. On the same graph we also plot the to- 
tal stellar radii of the two protostars, calculated from our 
stellar evolution code given the measured accretion rates of 
the sink particles. The binary orbit of the two sink particles 
is a very close one, and their separation is frequently com- 
parable to their total radius. In our original simulations, the 
sinks were treated as point masses once formed. This is still 
a good approximation for the majority of the sinks in our 
simulations which do not have such close encounters. How- 
ever, for extreme cases such as that shown in Figure (5] this 
assumption clearly breaks down. 
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Table 2. The number of sink pairs where the sinks pass within 
each other's tidal radius, compared to the total number of sink 
particles formed in each halo simulated in Sll. 



Halo 1 2 3 4 5 



Within Tidal Radius 


3/11 


0/8 


2/6 


1/6 


5/23 


Radii Touching 


2/11 


0/8 


2/6 


1/6 


1/23 



The tidal radius for a fluid satellite was calculated by 
Roche in 1847 to be 

n^2A4R^^ezy\ (6) 

where r* is the tidal radius, Rp the primary radius, pp the 
primary density and ps is the density of the secondary. The 
protostars considered here are well approximated as fluids as 
the outer extended layers will be easily separated from their 
parent protostar. Protostars which have encounters are usu- 
ally formed close to each other and have similar masses, radii 
and accretion histories. Consequently their internal densi- 
ties are also likely to be similar. Therefore we should expect 
tidal disruption of protostars which come within a distance 
of about 2.5 times the primary radius. To estimate the num- 
ber of protostars in each halo likely to be affected by tidal 
effects we calculate the number of protostars which come 
within a distance of 3 primary radii of each other (we round 
up to the nearest integer to account for the effects of the sec- 
ondary mass and oblateness). Since we do not have the full 
stellar evolution calculation for each star in our halo we use 
the radius from the semi-analytic model used in Sll which 
we introduced in Section |3] 

Table [2] shows the results. In all of the halos, bar one, 
there will be tidal disruption or mass transfer between a 
few protostars. A comparison to local star formation which 
might be of use to understa nd the physics of su ch an en- 
counter is the recent work of lLaioie fc Silld (|201lh in which 
mass transfer is seen between an approximately equal mass 
binary on an elliptical orbit after the system has reached 
periastron. We also calculate in how many cases the stel- 
lar radii go on to touch (see Table ^ , and see that this is 
also common. We can therefore conclude that mergers will 
contribute to the growth of some Population III stars, al- 
though they will not necessarily be the dominant mode of 
protostellar growth. 

Protostars that form a binary without merging may un- 
dergo mass transfer via Roche lobe overflow. This will in- 
fluence their subsequent evolu tion and orbital separation. 
iKrumholz fc ThompsonI l|2007l ) examine such an evolution 
in the context of present-day massive star formation, where 
the protostars have extended radii like those seen here due 
to high accretion rates. They show that in many cases the 
binary stellar mass ratio approaches unity and the orbital 
separation shrinks owing to the rapid mass exchange. The 
final outcome is a tight (< lAU) massive binary consist- 
ing of two almost identical stars. Such 'massi ve twins' are a 
comm on occurrence in the local universe (e.g. lHilditch et al] 
and if this were also true in primordial star forma- 
tion, they would be strong candidates for the progenitors 
of gamma-ray bursts (GRBs) in t he ea rly universe (e.g. 
iBromm fc Loebll2006l ). iFrver et all l|l999l ) have shown that 




L_, , , ^ I , , ^ , I , ^ , , I ^ , , , 

5 10 15 20 

Stellar Mass [M,„|] 



Figure 6. Effects of cold (low-entropy) mass accretion on the 
evolution of the stellar radius with the same variable accretion 
history as in Figure 2. The blue and green lines represent the cases 
where cold accretion is adopted after the stellar mass exceeds 
5 Mq and 2.5 Mq. The red line represents the fiducial case, where 
hot accretion is assumed throughout the evolution. 



evolved binary systems consisting of compact objects, such 
as neutron stars or black holes, can lead to GRBs from ac- 
cretion (Wooslcy 1993), or through mergers. Such systems 
could also form X-ray binaries during Roche lobe overflow 
from a post-main-sequence secondary on to a compact pri- 
mary, and in the ca se of mergers ma y even be a source of 
gravitational waves (|Frver et al.ll 19991 ). 

While not all protostars will merge or undergo mass 
transfer, their extended fluffy structure introduces addi- 
tional degrees of freedom to the stellar system during dy- 
namical interactions. This means that energy can be dis- 
sipated during close encounters by tidal distortion of the 
surfaces of the protostars. As a result, the post-encounter 
velocities of protostars that undergo close encounters may 
be smaller than the values computed in current simulations 
in which the sink particle is approximated as a solid body. 
Future studies of fragmentation in primordial halos should 
attempt to include these effects to obtain the most accu- 
rate description of the growth of Population III protostars 
and their flnal masses. Recentlv lCreif et al.l (|2012l ) made the 
flrst step in this direction in a model of the first 10 yr of a 
protostellar disks evolution. 



6 EFFECTS OF COLD MASS ACCRETION 

So far we have focused on the evolution of protostars under- 
going thermally efficient mass accretion ("hot accretion"). 
As discussed in Section 3, hot accretion is the most ap- 
propriate model for the rapid mass accretion expected in 
Population III star formation. However, for a better un- 
derstanding of the protostellar evolution, it is useful to 
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examine the effect of reducing the thermal efficiency of 
the accretion. E xperiments of this kind were presented in 
iHosokawa et al.l (2011a), where we examined the potential 
effects of different thermal efficiencies by switching the outer 
boundary conditions at some point during the evolution. 
Figure [6] shows the evolution of the stellar radius in such 
cases, taking the same variable accretion history as in Fig- 
ure 2. 

The evolution depends critically on the protostellar 
mass, A'h, at the point at which the boundary conditions 
are switched. With the cold accretion for M, > 5 Mq, the 
maximum stellar radius is almost the same as in the fiducial 
hot accretion case. When the boundary condition is switched 
at M« = 2.5 Mq, the protostar initially contracts and then 
re-expands to a radius of ~ 100 Rq at a mass of around 
10 Mq. This indicates that the cold mass accretion reduces 
the maximum stellar radius. The explanation for this phe- 
nomena comes from the stellar radius being at its maximum 
just prior to the KH contraction stage. In this transitional 
phase, heat accumulated in the stellar interior is transported 
outward and escapes from the star. A gas layer near the stel- 
lar surface will receive a portion o f this transported en tropy, 
which makes the star infiate (e.g.. IStahler et al]|l986l ). This 
effect would be alleviated if the mass accretion was thermally 
inefficient and the average entropy in the stellar interior was 
reduced. 



7 CONCLUSIONS 

We have combined the accretion rates measured from Pop- 
ulation III protostars in hydrodynamic simulations with a 
stellar evolution code to determine the structure of the pro- 
tostars. The measured accretion rates had a large degree of 
variability and were particularly high during the adiabatic 
stage of the protostars' evolution. This increased the stellar 
radius of the protostars, causing them to become puffy ex- 
tended objects. In the early stages of the protostar's life this 
intermittently in creased the a c cretio n rate of the protostar, 
but as shown in ISmith et al.l (|201ll ) this has only a small 
effect on the dynamic evolution of the halo. In some cases 
the separation of the protostars during their dynamic evo- 
lution was sufficiently close that mass transfer and mergers 
may occur between protostars. Mass transfer between mas- 
sive protostellar systems could lead to approximately equal- 
mass massive binaries that are a possible precursor of GRBs. 
The extended size and relatively large tidal radii of Popu- 
lation III protostars implies that interactions between them 
are likely to dissipate more energy than if they were pure 
point masses. 
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